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Migration cn profondeur 
avant addition de.donnees 



(54) Perfection nements aux precedes d'inversion tomographique d'evenements pointes sur les 
donnees sismiques migrees 

(57) Procede d'inversion d'evenements pointe sur la 
sismique migree en profondeur, caracterise en ce qu'on 
determine par migration profondeur avant addition une 
collection de traces se reflechissant a I'aplomb d'un 
point de surface donne et on met en oeuvre sur cette 
collection de traces un traitement qui selectionne une 
parametrisation du champ de vitesse qui optimise I'aii- 
gnement de points qui sont migres a partir de ladite col- 
lection de traces. 

Pointage d'une profondeur et d'un 
pendage de reflecteur au droit 
d'un point de surface (x, y). 



Demigration : 
S, R,t/c* , Vt/« independants du modele 
de vitesse. 
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Inversion par migration en profondeur . 
Determination du champ de v'rtesses pour 
leque! les points remigres sont alignes. 



FIG.1 



Printed by Jouve, 75001 PARIS (FR) 



BNSDOClD:<EP 1 239304 A1 I > 



1 



EP 1 239 304 A1 



2 



Description 

[0001] La presente invention est relative aux prece- 
des conversion tomographique d'evenements pointes 
sur des traces sismiques. 

[0002] Pour reconstruire des images du sous-sol, les 
geologues ou geophysiciens utilisent classiquement 
des emetteurs acoustiques places par exemple en sur- 
face. Ces emetteurs emettent des ondes qui se propa- 
gent dans le sous-sol et se reflechissent sur les surfaces 
des differentes couches de celui-ci (reflecteurs). Les on- 
des acoustiques reflechies vers la surface sont enregis- 
trees en fonction du temps par des recepteurs. Les si- 
gnaux enregistres par les recepteurs sont appeles tra- 
ces sismiques. 

[0003] On salt classiquement pointer sur de telles tra- 
ces sismiques des portions qui correspondent a la re- 
flexion d'impulsions acoustiques emises en surface et 
qui correspondent a des reflecteurs d'interet, ainsi que 
determiner les temps de trajet qui correspondent a ces 
reflexions. 

[0004] Les techniques d'inversion tomographique 
consistent a determiner, en fonction des traces sismi- 
ques acquises et d'une selection d'evenements, des 
moderations des champs de Vitesse dans le sous-sol. 
[0005] De nombreuses techniques d'inversion sont 
deja connues. 

[0006] II est en particulier classiquement connu d'in- 
verser les temps d'arrivee des ondes reflechies et no- 
tamment de determiner le champ de vitesse dans le 
sous-sol en minimisant I'ecart entre les temps d'arrivee 
observes et les temps d'arrivee modelises, ces derniers 
etant calcules par trace de rais, dans le modeie d'inver- 
sion, entre les emetteurs et les sources. 
[0007] L'inversion des temps de trajets s'avere diffici- 
le a mettre en oeuvre dans le cas de structures geolo- 
giques complexes, du fait notamment d'un manque de 
discrimination des arrivees multiples. 
[0008] D'autres techniques d'inversion tomographi- 
que sont des techniques qui mettent en oeuvre des mi- 
grations en profondeur. 

[0009] Une migration en profondeur avant addition 
des donnees sismiques consiste a determiner pour cha- 
que position x,y de surface une collection de traces mi- 
grees portant des evenements qui decrivent le sous-sol 
a I'aplomb de ladite position de surface (x, y). Une telle 
collection de traces est generalement rangee par clas- 
ses de distances source-recepteur (distances egale- 
ment appelees distances de deport ou d'offset) et ega- 
lement par classes d'angles speculates croissants, et 
plus rarement par classes d'orientation du segment 
source- recepteur. 

[001 0] Lorsque le modeie de vitesse utilise ref lete cor- 
rectement la realite, la profondeur associee a un eve- 
nement se reflechissant a I'aplomb d'une position de 
surface est sensiblement constante quelle que soit la 
distance de deport ou I'orientation de la trace conside- 
ree (ou encore Tangle speculate). 



[0011] Des procedes d'inversion ' tomographique en 
profondeur utilisent cette caracteristique pour definir le 
critere d'inversion a minimiser. 
[0012] Notamment, il a ete propose dans : 

5 

- "Velocity analysis by iterative profile migration" - 
Kama! Al-Yahya - Geophysics, Vol. 54, n°6, 1989, 
p. 718-729 

un procede d'inversion dans lequel on determine le 
10 modeie de vitesse du sous-sol en minimisant I'ecart 
a I'horizontalite de la courbe profondeur/distance 
d'une collection de traces migrees. 

[001 3] La mise en oeuvre de cette technique suppose 
is que plusleurs migrations en profondeur soient effec- 
tueesi ce qui est particulierementfastidieux a mettre en 
oeuvre lorsque plusieurs iterations sont necessaires 
pour converger. 

[0014] Un but del'invention est de proposer une tech- 
no nique d'inversion qui reduit sensiblement le nombre de 
migrations profondeur a effectuer pour construire un 
modeie optimum et ne necessite pas une puissance de 
calcul importante. 

[001 5] Elle propose a cet effet un procede pour la mi- 
25 se a jour d'un modeie de vitesse d'un sous-sol, selon 
lequel : 



on met en oeuvre sur un ensemble de traces sismi- 
ques acquises au droit dudit sous-sol une migration 
profondeur avant addition des donnees, cette mi- 
gration etant mise en oeuvre a I'aide d'au moins un 
modeie de vitesse de depart et permettant de de- 
terminer une ou plusieurs collection de traces qui 
decrivent chacune ie sous-sol a I'aplomb d'un point 
de surface, 

on pointe sur au moins une collection de traces ob- 
tenue a I'aide de cette migration un ou plusieurs 
evenements qui se reflechissent a I'aplomb du point 
de surface considere et on determine pour chaque 
evenement pointe une profondeur et un pendage 
de reflecteur et un pendage de reflecteur au droit 
dudit point de surface, 

on met eh oeuvre un trace de rais entre ce reflecteur 
et la surface pour determiner des collections de 
couples source acbustique-recepteur, ainsi que" 
donnees caracteristiques des temps de trajet et des 
gradients temps qui correspondent aux traces de 
rais associes a ces couples, 
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50 caracterise en ce qu'on met ensuite en oeuvre un 

traitement selon lequel on itere les etapes consistant a : 



migrer cinematiquement des donnees de temps de 
trajet et de gradient temps obtenues preced em- 
inent a I'aide d'une parametrisation du champ de 
vitesse, 

caracteriser I'alignement des points migres ainsi 
obtenus, 
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mettre a jour la parametrisation, 

le traitement selectionnant la parametrisation de champ 
de vitesse qui optimise I'alignement desdits points mi- 
gres. 

[0016] Les trois etapes precedentes peuvent etre re- 
petees jusqu'a I'obtention d'un alignement juge suffi- 
sant. 

[0017] On notera que le traitement de selection de la 
parametrisation du champ de vitesse ainsi propose ne 
necessite pas d'effectuersystematiquement une migra- 
tion profondeur des donnees sismiques (fort couteuse) 
apres chaque mise a jour des vitesses. 
[0018] On comprend par consequent que le procede 
propose par ^invention peut etre mis en oeuvre sans ne- 
cessiter des puissances de calcul importantes. 
[0019] Dans une mise en oeuvre avantageuse, la mi- 
gration en profondeur avant addition des donnees est 
mise en oeuvre pour une plurality de champs de vites- 
ses (CRP scan par exemple) et pour determiner la pro- 
fondeur et le pendage du reflecteur, on selectionne celui 
de ces champs de vitesses qui minimise I'ecart a I'hori- 
zontalite de la courbe profondeur/distance de deport ou 
de la courbe profondeur/angle speculate. 
[0020] Par ailleurs, la profondeur de reflecteur peut 
etre avantageusement determinee pour differents de- 
ports ou angles speculates. 

[0021] Selon un mode de mise en oeuvre prefere pos- 
sible, pour selectionner une parametrisation du champ 
de vitesse qui optimise I'alignement de points qui sont 
migres a partir des donnees de temps et de gradient 
temps, on -determine la distance entre les portions de 
reflecteur vues pour differents deports de la collection, 
la parametrisation de champ de vitesse selectionnee 
etant celle qui minimise cette distance. 
[0022] Selon un autre mode de mise en oeuvre, pour 
selectionner une parametrisation du champ de vitesse 
qui optimise I'alignement de points qui sont migres a 
partir des donnees de temps et de gradient temps, on 
determine les differences de profondeurs entre des 
points d'interception entre une verticale donnee et des 
portions de reflecteur qui sont paralleles au reflecteur 
mis a jour precedemment et qui passent par les nou- 
veaux points migres, la parametrisation de champ de 
vitesse selectionnee etant celle qui minimise s cette dif- 
ference. 

[0023] En variante encore, les pendages considered 
pour chaque deport d'une meme collection peuvent ne 
pas etre paralleles, auquel cas ils peuvent etre determi- 
nes a partir des rais. 

[0024] Egalement, on peut mettre en oeuvre les eta- 
pes suivantes pour selectionner une parametrisation du 
champ de vitesse, : 

on determine pour un point migre un reflecteur pas- 
sant par celui-ci, 

on trace lesrais entre ce point migre et la surface et 
on determine une erreurtemporelle qui est fonction 



de I'ecart entre le temps de trajet correspondent a 
au moins un de ces rais et le temps de trajet de la 
trace qui dans la collection correspond a la meme 
distance de deport ou au meme angle speculate, 
5 - on selectionne une parametrisation de champ de 
vitesse qui minimise cette erreur temporelle et op- 
timise Talignement de points qui sont migres a partir 
des donnees temps et de gradient temps. 

10 [0025] En variante ou en complement, on peut egale- 
ment mettre en oeuvre les etapes suivantes : 

on determine pour un point migr§ un reflecteur pas- 
sant par celui-ci, 
*5 - on trace les rais entre ce point migre et la surface et 
on determine une erreur temporelle qui est fonction 
de I'ecart entre le gradient de temps de trajet cor- 
respondent a au moins un de ces rais et le gradient 
de temps de trajet de la trace qui dans la collection 
20 correspond a la meme distance de deport ou au me- 
me angle speculate, 

on selectionne une parametrisation de champ de 
vitesse qui minimise cette erreurtemporelle et op- 
timise I'alignement de points qui sont migres a partir 
25 des donnees de temps et de gradient temps. 

[0026] D'autres caracteristiques et avantages de I'in- 
vention ressortiront encore de la description qui suit, la- 
quelle est purement illustrative et non limitative et doit 
30 etre lue en regard des dessins annexes sur lesquels : 

la figure 1 est un organigramme qui illustre un mode 
de mise en oeuvre possible pour I'invention ; 
la figure 2 est une representation schematique d'un 
35 sous-sol illustrant la mise en oeuvre de certaines 
etapes du traitement de la figure 1 ; 
la figure 3 est une representation schematique illus- 
trant d'autres etapes. 

40 [0027] Le traitement illustre sur la figure 1 comporte 
une premiere etape 1 qui consiste a migrer en profon- 
deur les donnees sismiques pour determiner pour un 
point de surface de coordonnees X, Y donnees les eve- 
nements qui sur une collection de traces sismiques cor- 

'45- ^respondent a une reflexion a I'aplomb de ce point. 

[0028] Cette migration est determinee avec un mode- 
le initial de vitesse de sous sol. Ce modele peut etre 
constitue en plusieurs couches. 

[0029] Pour chacune des couches de ce modele ini- 
so tial, le champ de vitesse est par exemple du type 



V(x,y,z)= (V0(x,y)+k(x,y) *z)*An(£,8,0phase) 

55 ou: 

x, y et z sont les coordonnees du point du sous-sol 
que Ton considere, z etant la profondeur ; 
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V(x, y, z) est la vitesse acoustique au point de coor- 
donnees x, y, z du sous-sol ; 
V0(x, y) et k(x, y) sont deux fonctions qui dependent 
des coordonnees x et y 

An(e(x, y), 5(x, y), 6) etant un terme decrivant I'ani- 5 
sotropie et pouvant dependre de x et y, ou £,5 sont 
par exemple les parametres d'anisotropie faible de- 
crits par Thomsen (Geophysics 1986), et ou 9 est 
i'angle de phase decrivant la direction de propaga- 
tion des ondes. 

[0030] Le traitement de migration profondeur avant 
addition peut egalement avantageusement etre mis en 
oeuvre pour une pluralite de modeles de vitesse corres- 
pondent a differentes perturbations possibles par rap- 
port au modele de vitesse initial: 
[0031] Ces perturbations sont par exemple de type li- 
neaire etfonction d'un coefficient de perturbation auquel 
on donne differentes valeurs pour disposer de differents 
modeles voisins du modele de champ de vitesse initial. 
La perturbation peut s'appliquer soit a une seule cou- 
che, soit a I'ensemble du modele. 
[0032] Pour reprendre I'exemple de modele de Vites- 
se initial donne precedemment, un modele de vitesse 
perturb6 pourra etre : 

V a (x,y,zHaVO(x,yH K(x } y)*z)*An 

dans le cas d'une perturbation appliquee uniquement a 
une couche ou 

V a (x,y,z)=a*(VO(x,y)+ K(x,y)*z)*An 

dans le cas d'une perturbation appliquee uniquement a 
I'ensemble du modele, 

ou a est le parametre de perturbation et a une valeur 
qui varie par exemple de 0,8 a 1 ,2. 
[0033] Cette technique de balayage de plusieurs mo- 
deles de vitesse au moyen d'un coefficient lineaire est 
appelee "CRP-scan" et a par exemple ete decrite dans 
le cas de migration 3D dans I'article suivant : 

Audebert, Diet, Zhang, "GRP* scans for 3D pre- 
stack Depth Migration : A powerful combination of 
CRP-gathers and velocity scan", SEG, 1996, Ex- 
panded Abstracts - p.515-518. 

[0034] Apres migration, on determine dans une 
deuxieme etape (etape 2 sur la figure 1) la profondeur 
et le pendage du reflecteur que Ton cherche a mettre a 
jour a I'aplomb du point de surface de coordonnees x, 
y considered 

[0035] Dans le cas ou plusieurs modeles de vitesses 
ont ete balayes, on pointe celui de ces modeles qui cor- 
respond a I'ecart a I'horizontalite de la courbe profon- 
deur/distance d'offset le plus faible, la profondeur Zee et 



le pendage Da que Ton determine pour le reflecteur 
etant alors fonction du modele ainsi pointe. 
[0036] Une fois ce pendage et cette profondeur de- 
terminee, on met en oeuvre avec le modele de vitesse 
considere (celui du coefficient lineaire a qui correspond 
a I'ecart a I'horizontalite le plus faible) un traitement in- 
verse d'un traitement de migration. 
[0037] Ce traitement - appele traitement de demigra- 
tion (figure 2) - consiste en des traces de rais pourcha- 
que deport et en fonction de Porientation source-recep- 
teur a partir d'un reflecteur (reflecteur Ref sur la figure) 
qui correspond rait a la profondeur (fonction du deport) 
et au pendage que Ton a determine a I'aplomb du point 
de surface de coordonnees x, y. 

[0038] Les donnees obtenues par ce trace de rais in- 
verse du reflecteur vers la surface ont I'avantage d'etre 
totalement independantes du modele de vitesse consi- 
dere. Ces donnees sont notamment les positions a la 
surface des couples source-recepteur (S, R) correspon- 
dant a une reflexion sur le reflecteur Ref, ainsi que les 
temps de trajets et les gradients temps qui correspon- 
dent aux rais calcules, ces temps de trajet et gradients 
temps etant independants du modele de vitesse consi- 
dere lorsqu'ils sont rap portes au coefficient a dans le- 
quel le traitement de demigration a ete mis en oeuvre 
(t/a et Vt/a sur la figure 2 dans le cas d'une perturbation 
appliquee au modele). 

[0039] On pourra a cet egard avantageusement se 
rapporter a I'article : Audebert, Diet, Guillaume, Jones, 
Zhang - CRP-scans : "3D PreSDM Velocity Analysis via 
Zero Offset Tomographic Inversion" - 1997 - Expanded 
Abstracts - Soc. Expl. Geophys. -pp. 1805-1808. 
[0040] Dans une etape suivante (etape 4 sur la figure 
1), on met en oeuvre avec les donnees ainsi obtenues 
un traitement consistant, pour une distance d'offset h 
choisie, a migrer par trace de rais pour determiner le 
point de reflexion du couple source/recepteur S, R qui 
correspond a cette distance d'offset h et qui honore les 
temps de trajet et les gradients temps pn§alablement 
constitues. 

[0041] Cette migration cinematique est faite soit dans 
le modele de migration, soit apres avoir introduit une 
perturbation sur le modele de vitesse utilise dans Petape 
3 de demigration. 

[0042] Comme I'illustre la figure 3, cette- nouve lie mi- 
gration permet de determiner un nouveau point migre 
Pmi, ce point etant defini par sa profondeur Zi, ainsi que 
par les coordonnees Xi, Yi qui correspondent aux coor- 
donnees du point de surface a I'aplomb duquel ce point 
Pmi se trouve. 

[0043] Cette migration est mise en oeuvre pour diffe- 
rentes distances d'offset - ou differents angles speculai- 
res, ce qui est equivalent - de la collection de traces ob- 
tenue a Tissue de I'etape 3. 

[0044] Comme I'illustre la figure 3, les differents 
points migres (en I'occurrence Pmi, Pm'i) obtenus pour 
ces differentes distances d'offset (en Poccurrence h, h 1 ) 
- ou ces differents angles speculates - sorft disperses 
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lateralement et ne sont pas concentres sur le point Pa. 
[0045] Le proc6de propose tire toutefois partie de ce 
que lorsque le modele de Vitesse que Ton utilise reflete 
correctement la r6alite, ces differents points migres doi- 
vent etre alignes sur une meme portion de reflecteur. 
[0046] Or, en faisant Thypothese que le reflecteur i 
que Ton cherche a mettre a jour peut se comparer au 
reflecteur j qui correspond a la couche mise a jour pre- 
cedemment (ce qui est d'autant plus proche de la reality 
que la courbure residuelle, est faible), il est possible de 
caracteriser le fait que les points remigres (Pmi, Pmi' 
sur la figure 3) sont plus ou moins alignes sur une meme 
portion de reflecteur. 

[0047] En faisant cette hypothese, on peut par exem- 
ple determiner le champ de vitesse que Ton recherche 
en minimisant la distance entre les portions de reflecteur 
vues par differentes traces de la collection migree dans 
I'etape 3 (distance D sur la figure 3), ces portions de 
reflecteurs etant supposees paralleies ou quasi paralle- 
ls au reflecteur j. 

[0048] Un autre critere encore peut consister a mini- 
miser les ecarts entre les profondeurs des points d'in- 
tercepts entre une verticale (par exemple celle a 
I'aplomb du barycentre des lieux migres) et les portions 
de reflecteur plus ou moins paralleies au reflecteur j mis 
a jour precedemment et passant par les nouveaux 
points migres (ecart entre les coordonnees Zi, Z'i sur la 
figure 3). - 

[0049] Dans I'un et I'autre cas, les minimisations sont 
mises en oeuvre en modifiant par iterations successives 
les parametrisations du modele de vitesse. 
[0050] A titre d'exemples, les perturbations de para- 
metrisation utilisees successivement peuvent etre de- 
terminees en fonction des traces de rais calcules lors 
de chaque nouvelle migration, pour les distances d'bff- 
set ou les angles speculates consideres, de la fagon 
qui est proposee dans : 

Liu et al., 1 994 - "Velocity analysis by perturbation" 
- Sixty fourth Annual International Meeting, Soc. Ex- 
pl. Geophys., expended abstracts, 1991-1994; 
Liu etal., 1995 -"Migration velocity analysis : theory 
and iterative algorithm", Geophysics - 60-142-1 53 ; 
Wang et al. - "Macro velocity model estimation 
through model based globally optimised residual 
curvature analysis" - Expended abstracts - Soc. Ex- 
pl. Geophys., 1596-1599 (1998), 
les techniques decrites dans ces articles utilisant le 
fait que les champs de vitesse v et les profondeurs 
z des points de reflexion des rais verifient : 

3z/8v= A3t/8v 

ou t est le temps de trajet du rai et A un para- 
metre qui est fonction du pendage du sous-sol et 
est calculable notamment a partir du rai. 



[0051] Parailleurs, letraitementde selection de lapa- 
rametrisation de champ de vitesse peut prendre en con- 
sideration d'autres criteres en plus de I'optimisation de 
I'alignement de points qui sont migres a partir de ladite 
5 collection de traces. 

[0052] Notamment, il est egalement possible de de- 
terminer, par exemple pour la distance de deport qui cor- 
respond au point migre Pmi le plus proche du reflecteur 
Ref, la portion de reflecteur qui correspond a ce point 
■ Pmi et qui est parallele a la couche J mise a jour. Puis, 
pour chacune des traces de la collection obtenue a Tis- 
sue de I'etape 3, on trace le rayon qui se reflechit sur la 
portion de ce reflexion No. Les rais ainsi obtenus ont 
des temps de trajet et des gradients de temps differents 
de ceux obtenus a Tissue de I'etape 3. Les ecarts entre 
les temps ou gradients de temps obtenus a Tissue de 
P£tape 3 et ceux obtenus par ces traces de rais consti- 
tuent des attributs d'erreurs qui permettent egalement 
de caracteriser la convergence du champ de vitesse 
choisi par rapport au champ de vitesse reel. 



Revendications 

25 1. Procedepour la mise a jour d'un modele de vitesse 
d'un sous-sol, selon lequel : 

on met en oeuvre sur un ensemble de traces 
sismiques acquises au droit dudit sous-sol une 



30 migration profondeur avant addition des don- 

nees, cette migration etant mise en oeuvre a 
Taide d'au moins un modele de vitesse de de- 
part et permettant de determiner une ou plu- 
sieurs collection de traces qui decriventchacu- 

35 ne le sous-sol a I'aplomb d'un point de surface, 

on pointe sur au moins une collection de traces 
obtenue a Taide de cette migration un ou plu- 
sieurs evenements qui se reflechissent a 
I'aplomb du point de surface considere et on 
determine pour chaque evenement pointe une 
profondeur et un pendage de reflecteur et un 
pendage de reflecteur au droit dudit point de 
surface, 

on met en oeuvre un trace de rais entre ce re- 
45 " ' " '"^flecteur et la surface pour determiner des col- 
lections de couples source acoustique-recep- 
teur, ainsi que donnees caracteristiques des 
temps de trajet et des gradients temps qui cor- 
respondent aux traces de rais associes a ces 
50 couples, 

caracterise en ce qu'on met ensuite en 
oeuvre un traitement selon lequel on itere les eta- 
pes consistant a : 

55 

migrer cinematiquement des donnees de 
temps de trajet et de gradient temps obtenues 
precedemment a Taide d'une parametrisation 
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du champ de vitesse, - 

- caracteriser I'alignement des points migres ain- 
si obtenus, 

- mettre a jour la parametrisation, 

le traitement selectionnant la parametrisation de 
champ de vitesse qui optimise I'alignement desdits 
points migres. 

2. Precede selon la revendication 1 , caracterise en 
. ce que la migration en profondeur avant addition 
des donnees sismiques est mise en oeuvre pour 
une pluralite de champs de vitesses et en ce que 
pour determiner la profondeur et le pendage du re- 
f lecteur, on selectionne celui de ces champs de vi- 
tesses qui minimise I'ecart a I'horizontalite de la 
courbe profondeur/distance de deport ou de la 
courbe profondeur/angle speculaire. 

3. Procede selon I'une des revendications preceden- 20 
tes, caracterise en ce que la profondeur de reflec- 
teur est determinee pour differents deports ou an- 
gles speculates. 
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5. 



Precede selon la revendication 1, caracterise en 25 
ce que pour selectionner une parametrisation du 
champ de vitesse qui optimise I'alignement de 
points qui sont migres a partir des donnees sismi- 
ques de temps et de gradient temps, on determine 
la distance entre les portions de reflecteur vues 30 
pour differents deports, la parametrisation de 
champ de vitesse selectionnee etant celle qui mini- 
mise cette distance. 

Procede selon la revendication 1 , caracterise en 35 
ce que pour selectionner une parametrisation du 
champ de vitesse qui optimise I'alignement de 
points qui sont migres a partir des donnees de 
temps et de gradient temps, on determine les diffe- 
rences de profondeurs entre des points d'intercep- 40 
tion entre une verticale donnee et des portions cfe 
reflecteur qui sont paralleles au reflecteur mis a jour 
precedemment et qui passent par les nouveaux 
points migres, la parametrisation de champ de vi- 
tesse selectionnee-etant celle qui minimise cette 45 
difference.. 



6. Procede selon Tune des revendications preceden- 
tes, caracterise en ce que pour selectionner une 
parametrisation du champ de vitesse, : 50 

- on determine pour un point migre un reflecteur 
local passant par celui-ci, 

• on trace les autres rais de la collection entre ce 
reflecteur local et la surface et 55 

• on determine une erreur temporeUe qui est 
fonction de I'ecart entre le temps de trajet cor- 
respondant a au moins un de ces rais et le 



temps de trajet de la trace qui dans la collection 
correspond a la meme distance de deport ou 
au meme angle speculaire, 
- on selectionne une parametrisation de champ 
de vitesse qui minimise cette erreur temporelle 
et optimise ralignement de points qui sont mi- 
gres a partir des donnees sismiques temps et 
de gradient temps. 

Procede selon J'une des revendications preceden- 
tes, caracterise en ce que pour selectionner une 
parametrisation du champ de vitesse, : 

- on determine pour un point migre un reflecteur 
passant par celui-ci, 

- on trace les autres rais dela collection entre ce 
reflecteur et la surface et 

- on determine une erreur temporelle qui est 
fonction de I'ecart entre le gradient de temps 
de trajet correspondant a au moins un de ces 
rais et le gradient de temps de trajet de la trace 
qui dans la collection correspond a la meme 
distance de deport ou au meme angle specu- 
laire, 

- on selectionne une parametrisation de champ 
de vitesse qui minimise cette erreur temporelle 
et optimise I'alignement de points qui sont mi- 
gres a partir des donnees sismiques de temps 
et de gradient temps. 
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Migration en profondeur 
avant addition de donnees 



Pointage d'une profondeur et d'un 
pendage de reflecteur au droit 
d'un point de surface (x, y). 








Demigration : 
S, R, t/o( , Vt/« independants du modele 
de vitesse. 





Inversion par migration en profondeur . 
Determination du champ de vitesses pour 
lequel les points remigres sont alignes. 



FIG.1 
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Demigration 
modele o( 




FIG. 2 
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